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Abstract 

The binding energies and matter distributions for the 3- body sys- 
tems hke (p- meson +2A^, 2(f) + N and 4- body system hke (j) + 3n 
are calculated. For the 3- particle systems two- dimensional Faddeev 
equations in the differential form are used. For the 4- body system 
+ 3n the folding model is applied. 



1 Introduction 

As it has recently been intensively discussed [1-7], there are indications of 
strong attraction of mesons with one strange quark K~ (K) to few- nucleon 
nuclei. Along this line, it is interesting to look at the interaction of a meson 
with two strange quarks like the 0-meson with light nuclei. Already exist- 
ing theoretical investigations of the 0-meson interaction show rather strong 
attraction between a 0-meson and a nucleon. Indeed, the calculation of the 
(f) — N interaction within the quark model [S| and on the basis of a totally 
different phenomenological model [9j based on the dominant role of the ss 
configuration in the 0-meson structure predicts considerable (p — N attraction 
with a binding energy of about 9 MeV for the (pN system. 

Strong attraction is not very surprising in view of the physical arguments 
that the strong K~N attraction originates from the influence of subthreshold 
resonances A 1405 and S 1335 . 

Indeed, let us compare the mass of the state (p + N with the masses of two 
subthreshold states K + A 1405 and K + S1385. It turned out that distances 
of the above subthreshold states from the (f) + N threshold are of the same 
order of magnitude as in the K~N case, which means that as in K~N one 
can expect strong influence of A 1405 and S 1385 and strong attraction also in 
the (f)N system. 

Bearing in mind this sort of strong attraction in the (pN system, it is 
interesting to consider the possibility of bound states of a 0-meson with a 
few nucleons, in particular, with two neutrons or two protons. This is in fact 
a question concerning the existence of new nuclear clusters which do not exist 
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Figure 1: The 00 local potential 

without a (f)- meson. In what follows, we will calculate the binding energies 
of the three-body systems (pnn, (pnp and 0pp. 

2 Input 

It is reasonable to consider a 3- particle system of the type 00iV in the 
same theoretical frameworks. To describe this system, one needs apart from 
the (j) — N potential also the — potential. One of the simplest ways to 
construct this potential is to take it in the form of a sum of two Yukawa 
potentials (attractive and repulsive): 

V,,{r) ^V,^-V,^ (1) 

where 

Vi = lOOOMeV-fm, /xi = 2.5 fm-\ 

V2 = 1250 MeV-fm, /is = 3.0 fm"^ . 

The parameters of this potential were fixed by the position and width of 
the f2(2010) resonance which has only one decay channel into two 0- mesons. 
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Figure 2: The location of the pole of the absolute value of the — scattering 
amplitude in the complex energy plain. 

As in [9j , a Yukawa type potential is chosen for the (p — N interaction 

VMr) = -«e-^7r (2) 

with a = 1.25 and fi = 600 MeV. This potential is rather deep and narrow 
and supports binding in the (pN system with binding energies Ei/^n = —9.47 
and E^p = -9.40 MeV. 

For the np triplet s-wave interaction the potential MT III was used. 
Our singlet s-wave interaction is based on the potential MT I [TU] with a 
slight modification of having now a parameter = 2.617. This value is 
chosen in order to reproduce the experimental value of the nra-scattering 
length ann = —18.5 fm [11]. One can see that in the three-body systems 
(f)NN there are two scales of distance related to the different ranges of the 
N—N and (p — N interactions. This may produce a delicate interplay between 
a narrow attraction area of the cp — N interaction and repulsive parts of the 
MT-potentials, as it was emphasized in [12]. Apart from that, different 
ranges of the interaction can provide the cluster formation in the systems 
under consideration. 

Let us start to describe three particle systems. 
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3 Calculations 

Our calculations are based on Faddeev equations in a differential form 
written down for the 3-body systems (pNN and (jxpN. 

First, the Faddeev components of the wave function are expanded into 
partial waves: 

^ U^{Va,Ca)Y,^,°{fj^,U (3) 



where ?7„ = \f]a\, = Va = Va/\f]a\, L = ^a/ \^a\, are the 

bispherical harmonics. Jacobi coordinates f]a,^a were used, and only the 
lowest partial wave is taken into account. 
The Jacobi coordinates are as usual: 

Vi - r , = ^ (4) 

—* 

— — - = -r- 5 

rrii + rrij Oa 

where rj, denote the radius- vector and the mass of particle i, 



mirrii , Imdmi + m,) 

' ha = \ ^ ^,M = mi + m2 + m3 



V {mi + mj)M ' V 

and indices a take on the following values: a = 3 for {ij)k = (12)3, a = 1 
for {ij)k = (23)1, a = 2 for {ij)k = (31)2. 

Since there are two identical particles in the system (we take uin = rUn 
for (pnp system), the following two coupled-differential Faddeev equations 
survive: 

D + V,{^)- e] U,{p,v) = 



D + V, (^) - E ] U,{p, y.) 



(6) 
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{U3 = U, 



where polar coordinates p = ^rj^, + taiKpa = ^a/Va are introduced 
and 



D 



Vi — Vnn , V2 — Vtf)N , 

( Id 1 92 



2M y P dp p"^ dif"^ J 

C+ = Min{|v? + 7a'a| , Vr - (V9 + Ja'a)} , C- = \(p - 7a', 



7ij = arcsm Sij , Sij 



uikM 



V {rrii + ruk) (rrij + rUk) 
{ijk = 123, 231, 312) 

indices correspond to 1 for 0-meson, 2 and 3 for nucleons. 

The two-dimensional system of Faddeev equations has been solved by 
discretization of variables: hyperradius p and hyperangle if with N and M 
mesh points, respectively. Stable results for three digits of binding energies 
were reached at N = 110, M = 210 and L{p variable cutoff) = 9 fm. 

As a result, the binding energy of the system (pnn with value -E^^^ = 
—21.8 MeV and value E^np = —37.9 MeV for the binding of the (pnp system 
with np pair in triplet state have been obtained. It should be noticed that 
for this binding energy in the (pnp system both the main 0-meson decay 
channels on i^-mesons are closed. Let us comment the last value of energy 
which appeared rather large. From naive reasons in the configuration (p + d 
one would expect binding of an order of 2 x E^^ + E ^ which is much smaller 
than the calculated value. However, due to the strong attraction in the (pN 
- subsystem {E^i^n ~ —9 MeV) one can expect that in the 3-particle (pnp 
system the configuration (p + d is rather suppressed. Hence, it follows that 
in the above system there is no strong cancellation between potential and 
kinetic energies of nucleons, like in deuteron, and a strong attractive triplet 

- potential [Vt ~ 100 MeV) shows its full value. 

The dependence of the (pnn binding energy on the parameter a of the 
(j) — N interaction is investigated. It is shown in Figure [3] that excited states 
appear in this system. 

As can be seen from these results, the binding in 3-particle systems like 
(f)NN is possible even for weaker (p — N attraction than the one of poten- 
tial (121). 
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Now let us turn to the (j)(j)N system. The indication of a resonance in the 
+ subsystem and its parameters have been found by drawing the Argan 
diagram and studying the position of the pole in the scattering amplitude. 

With the potential V^^ of the form ([1]), the binding energy of the 3-body 
system 00n turns out to be - 77 MeV. 

As in the case of the (p{np)tripiet state, due to the strong binding in the 
00iV system, the K-meson decay channels of the - meson are closed. 

Let us discuss now the 4-body system 0nnn. To make a preliminary 
estimate of its binding energy, we have used the folding model with the 
(0nn)-subsystem cluster. 

By averaging the interactions of the third neutron with the particles of the 
(07zn) system over the cluster wavefunction it is easy to obtain an effective 
potential which has the form shown in Figure |H Here r means the distance 
between the third neutron and the center of mass of the {(pnn) cluster. 

The solution of the Schroedinger equation with this potential shows that 
there are no bound states. However, the exact treatment of an analogous 4- 
body system (rjc+SN), performed by means of the AGS- equations [TB], shows 
that the folding model may strongly underestimate the exact calculation. 

We plan to perform such a calculation in the near future. 

Finally, let us discuss the matter distribution, shown in Figures [5] and [6l 
for the 0nn system. First of all, as one can see on both the pictures which 
show the profiles of two Faddeev components of the wave function, the main 
part of a matter is concentrated at small distances inside the volume of 
the size of 2 fm. In addition, one can see that the fine structure of different 
components of the wave function is nonhomogeneous. Indeed, the component 
U2 has an extremum at hyperangle ~ 71/ 2 which means that a 0- meson 
and one of neutrons are close to each other forming the substructure in the 
(pnn system, which was mentioned at the end of the section 2. In contrast 
to that, the component of the wave function Ui has two extrema without 
well- pronounced substructures. Having in mind the small size of the system 
and its neutrality, one would expect the propagation of an object like this 
far away from the region of production. 
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Figure 3: The dependence of the binding energy of the (pnn system on the 
parameter a of the 4> — N interaction. 
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Figure 4: The folding potential with the p-wave centrifugal barrier for the 
four-body system {(pnn) + n. 
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Figure 5: The partial amplitude Ui{p, ip) for the (pnn system. Extreme values 
of Ui{p, (f) are 6.772 x 10 (at p = 0.654, (p = 1.092) and -5.633 x 10 (at 
p — 1.881, (p = 0.695). Contours of different cross-sections of the Ui{p,Lp) 
surface are shown in the (p, ip) plain. Dotted, dashcd-dottcd and second 
dashed-dotted hues correspond to vahies 5.643 x 10 1.693 x 10 —5.121 x 
10^'^ on the Ui{p,ip) axis respectively. 




Figure 6: The partial amplitude U2{p,(p) for the (f)nn system. The extreme 
value of U2{p,ip) is -2.182 x lO'^ (at p = 1.554, ip = 1.323). Contours of 
different cross-sections of the U2{p,p>) surface are shown in the {p,ip) plain. 
Dashed-dotted, dotted and dashed lines correspond to values —1.983 x 10 
— 1.454 X 10"^ and —0.545 x 10 on the U2{p,(p) axis respectively. 
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